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SYNOPSIS 


The  Ritz  energy  method  le  used  to  obtain  solutions  for  an 
interior  panel  of  a plate  or  slab  which  Is  continuous  over  a rectang- 
ular grid  of  flexible  beams  supported  by  columns  at  their  Inter- 
sections* It  Is  assumed  that  the  plate  Is  uniformly  loaded  oyer  its 
whole  area,  and  that  parallel  beams  are  of  equal  stiffness  end  are 
spaced  regularly.  The  plate  deflection  is  represented  by  an  infinite 
series  of  polynomial  functions  derived  by  U*  J.  Duncan  for  beams 
clamped  at  both  ends;  these  functions  lead  to  relatively  accurate 
results  with  a reasonable  amount  of  maerlcal  work*  Deflections, 
moments  and  average  moments  for  ceses  covering  a vide  range  of  beam- 
plate  rigidities  and  width-length  ratios  are  presented  In  tabular 
and  graphical  form* 
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I,  UTECDDCTICH 

1.  Object  and  Scop*  of  Investigation 

Sene  of  tha  most  Important  apodal  problems  of  Elastlo  Theory 
aro  thoao  involving  a medium-thick  rectangular  plato  aupportad  along  ita 
boundarioa  and  bant  by  loada  applied  to  one  of  ita  faces.  In  many  cases, 
particularly  thoao  in  which  two  oppoaito  adgaa  of  tha  plate  are  simply 
supported , theoretically  cample to  solutions  of  the  problems  hare  been 
obtained.  However,  for  many  other  cases  exact  solutions  have  not  been 
found;  the  most  notable  example  is  the  case  of  a uniformly  loaded  plate 
fixed  against  displacement  and  rotation  at  its  edges.  In  studying  such 
problems  recourse  has  been  had  to  experimental  methods^  or  to  approx- 
imate numerical  methods  of  solution.  In  the  present  work,  an  approx- 
imate numerical  solution  is  obtained  for  the  technically  important  case 
of  a plate  supported  by  elastic  beams. 

The  problem  considered  is  that  of  an  interior  panel  of  a plate 
continuous  over  a rectangular  grid  of  flexible  beams  which  are  supported 
at  their  intersections  by  columns ^ It  is  assumed  that  parallel  beams 
are  of  equal  stiffhess  and  are  uniformly  spaced  and  that  the  system  has 
a large  number  of  panels  in  both  directions.  The  plate  is  considered  to 
be  acted  upon  by  a lateral  load  whioh  is  uniformly  distributed  over  its 
whole  area.  Qader  these  conditions,  the  bending  in  all  interior  panels 
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For  instance  for  the  clamped  plate  problem  see  references  (1)  and  (2) 
of  the  Bibliography. 


2 


nay  be  assumed  to  to  identical,  consequently  only  ona  interior  panel  need 
to  considered.  Apart  from  their  intrinsic  Importance 9 results  of  the 
analysis  for  this  case  of  uniform  loading  of  the  entire  plate  can  to  used 
to  obtain  moments 9 shears  and  deflections  for  several  other  loadings  by 
superposition  with  known  solutions. 

The  problem  described  has  been  solved  by  the  use  of  a method 
which  embodies  the  customary  assumptions  for  the  flexure  of  medium-thick 
plates.  In  brief  these  are: 

1.  The  plate  is  of  uniform  thickness  and  is  homogeneous  and 
l8otroplcv  at  least  so  far  as  strains  and  defamations 
parallel  to  the  bounding  surfaces  are  concerned. 

2.  Hooke’s  law  applies  to  the  horizontal  strains 9 which  are 
assumed  to  vary  linearly  through  the  thickness.  This 
implies  that  the  deformation  due  to  shear  is  negligible 
compared  with  that  due  to  the  direct  bending  stresses* 

3.  Compressive  stresses  between  planes  parallel  to  the  bound- 
ing surfaces  are  small  compared  to  the  flexural  stresses 
and  may  to  neglected. 

4.  The  resultant  of  the  direct  stresses  acting  on  any  cross- 
sooti  on  of  the  plate  is  assumed  to  to  a couple.  It 
can  to  shown  (3)*  that  this  is  equivalent  to  requiring 
that  the  deflections  of  the  plate  (measured  from  a develop- 
able surface)  be  small  in  comparison  with  the  plate  thick- 
ness. These  assumptions  follow  very  closely  those  which 
•re  customary  In  the  ordinary  theory  of  flexure  of  beams. 


* Bombers  in  parentheses 9 taxless  otherwise  identified 9 refer  to  the 
Bibliography. 


3* 


In  addition  to  these  limitation*  tha  problem  has  baan  simpli- 
flad  by  assuming  that  tha  widths  of  tha  beams  and  tha  cross-sectional 
dimensions  of  columns  ara  compared  vith  tha  paaal  dimensions*  Thia 

simplification  will  not  have  much  affact  on  tha  daflactlona  and  moments 
naar  tha  eantara  of  tha  panala  but,  in  sobs  oases,  will  land  to  excessive- 
ly  high  moments  naar  tha  columns* 

Tha  Boat  cannon  structure  similar  to  tha  plate-beam  ayatan 
da sc ri bad  is  tha  two-way  concrete  slab*  This  type  of  structure  differs 
from  tha  idealised  problem  in  that  reinforced  concrete  is  neither  homo- 
geneous nor  isotropic  nor  perfectly  elastic*  Furthermore  the  supporting 
beams  9 whan  of  concrete  , a re  usually  built  nonollthically  vith  tha  slab* 
Tha  affact  of  casting  tha  slab  oonolithically  vith  tha  supporting  beams 
can  be  approximately  accounted  for,  in  tha  usual  manner,  by  assuming 
that  part  of  tha  slab  is  effective  as  a flange  on  tha  beam*  In  general v 
tha  ordinary  isotropic  theory  of  flexure  of  plates,  together  vith  ex- 
perimental data,  has  provided  a satisfactory  basis  for  design  of  slab 
elements  of  present  d*y  reinforced  concrete  structures* 

Two-way  concrete  slabs  are  usually  designed  without  accounting 
for  the  flexibility  of  the  supporting  beams,  that  is,  the  beam  deflections 
are  usually  assumed  to  be  small  in  comparison  with  the  deflections  of 
the  central  portions  of  the  loaded  panels*  The  present  work  gives  the 
theoretical  distribution  of  momenta  in  an  interior  panel  for  supporting 
beams  of  any  stiffness*  The  solutions  include,  as  limiting  cases  of  xhe 
general  problem,  the  known  results  for  plates  fixed  against  rotation  and 
*4  thaiT*  rninrwrt*.  and  tnr  nlaiaa  amiuorted  bv  a reotamffular 

*•  * ~ " m r ~ — ' • a — A ■ ■ - - -a  « * w 

array  of  columns  without  connecting  beams  (flat  slabs). 
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2.  Aefcnovladgnant 

This  inraatigation  has  baaa  la  part  snpportad  by  tha  Mechanics 
Branch,  Qffiea  of  laral  Base  arch,  Dapartnant  of  tha  Ravy,  in  connection 
with  a ra  tsa arch  pro  gran  antitlad  "Runarical  and  Approadnata  Methods  of 
Strass  Analysis, " in  prograss  in  tha  Structural  itasaareh  Laboratory,  of 
tha  Dapartnant  of  Civil  Enginaarlng,  and  in  tha  Enginaaring  Exparinant 
Station  of  tha  Unirarsity  of  Illinois,  undar  Contract  B&ori  - 71,  Task 
ordar  6.  Tha  natarial  in  this  raport  has  haan  drawn  frcn  a dissartation 
of  tha  sana  titla  suhnittad  to  tha  Qr adnata  Facility  of  tha  thivarsity  of 
Illinois  by  Mr*  J*  G*  Sutharland,  January  1953* 
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IT.  METHOD  OP  SOLUTION 


3*  Notation 

Tha  following  nomenclature  is  used',  Any  consistant  sat  of 
unite  can  be  employed)  hove re r,  for  refaranoa , tha  unite  are  given  in 
tha  anginaaring  inch-pound-second  system. 


a = length  of  penol,  in  the  x-direction  (in) 

b = width  of  panel,  in  tha  y-direction  (in) 

c = b/a  = ratio  of  panel  width  to  panel  length 
(nondlaansional) 

m,n  = integers  ranging  from  one  to  infinity 

q = intensity  of  uniform  load  on  tha  plate  (lb/in  ) 
r,e  = integers;  particular  values  of  m and  n respectively 
v = lateral  deflection  of  the  plate  (in) 
x,y  = Cartesian  coordinates,  directed  along  the  beams  (in) 

z1»y1  = Cartesian  coordinates  with  their  origin  at  the 
center  of  the  panel  (in) 

D = flexural  rigidity  of  the  plate  (in  lb) 

E ,E.  = moduli  of  elasticity  of  the  beams  of  span  length 
a and  b,  respectively  (lb/in2) 

I#,I^  = second  momenta  of  area  of  the  croas-eections  of  the 
**  beams  of  span  length  a and  b,  respectively  (inA) 

= bending  moments  per  unit  length  of  sections  of  a 
y plate  perpendicular  to  x-  and  y-axes,  respectively; 
the  superscript  "°"  is  used  to  Indicate  that 
Poisson's  ratio  is  zero  (lb) 

= average  bending  moments  per  unit  length  on  various 
sections  of  the  plate;  the  superscript  "°"  is  used 
to  indicate  that  Poisson's  ratio  is  zero  (lb) 

M ,M.  = momenta  in  the  beams  of  span  length  a and  b,  res- 
pectively (in  lb) 
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= mcnent  coefficients  for  computing  average  moments 
for  any  value  of  Pole eon’s  ratio  (lb) 

= a set  of  ajiiiti  In  functions  (nondlmenslonal ) 

= strain  energy  due  to  bending  (in  lb) 

= total  potential  energy  of  one  panel  (In  lb) 

= potential  energy  of  the  load,  q,  far  one  panel  of 
the  system  (In  lb) 

= mre/2  = term  appearing  In  the  exact  solution  for  a 
plate  supported  by  columns  only  (nondlmenslonal) 

= coefficients  defining  the  deflected  shape  of  the 
plate  (nondlmenslonal) 

= E#I yaD  = ratio  of  beam  rigidity  to  plate  rigidity 
for  each  bean  parallel  to  the  x-axla  (nondlmenslonal) 

= R^I^/bD  = ratio  of  beam  rigidity  to  plate  rigidity 
for  each  bean  parallel  to  the  y-axls  (nondlmenslonal) 

= Poisson’s  ratio 

= x/a  = dimensionless  coordinate  In  the  x-directlon 
= y/b  = dimensionless  coordinate  In  the  y-directlon 


4*  B£sun6  of  the  Bits  Energy  Method  as  Applied  to  Plate  Equilibrium 
Problems 

By  suitably  restricting  the  deflected  shape  of  a continuous 

elastic  system,  It  can  be  reduced  to  a similar  but  simpler  system  having 

a finite  number  of  degrees  of  freedom.  This  Idea  was  first  used  by  Lord 

Rayleigh^  who  applied  It  in  approximate  calculations  of  the  fundamental 

frequencies  of  vibration  of  complicated  structures.  In  these  applications 

the  restrictions  art  equivalent  to  constraints  acting  on  the  DtniuWs. 

(5) 

V.  Ritz  later  developed  the  same  idea  in  a more  general  form;  he  dis- 
cussed in  detail  the  solution  of  the  problem  of  the  rectangular  plate 
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clamped  at  ita  boundaries*  the  solution  of  boundary-value  problems 
governed  by  Laplace 's  equation*  the  solution  of  second-order  differential 
equations  with  variable  coefficients,  and  the  solution  of  the  problem  of 
the  vibrating  string.  The  general  method  Is  known  as  the  Rsyleigh-Rltz 
energy  method;  In  applications  to  plate  equilibria*  problems  it  is 
usually  referred  to  as  the  Rltz  method.  The  Raylelgh-Bitz  method  has 
been  one  of  the  most  useful  approximate  numerical  methods  in  the  field 
of  structural  analysis;  it  has  yielded  satisfactory  solutions  for 
numerous  problems  In  equUltorltm*  buckling  and  vibration. 

The  analytical  equilibrium  problem  for  plates  is  that  of  find- 
ing a deflection  function  v(x*y)  which  satisfies  the  equation  V^w  = q/D 
within  the  boundary  of  the  plate  and  which  also  satisfies  certain  con- 
ditions on  the  boundary.  Rltz  observed  that  the  mathematical  problem 
is  the  sane  as  that  of  rendering  stationary  the  integral  expression  for 
the  total  potential  energy  of  the  system.  With  the  ordinary  theory  of 
flexure  for  plates  and  beams*  the  total  potential  energy*  V*  can  be  ex- 
pressed In  terms  of  the  deflection  function,  w*  and  Its  second  deriv- 
atives. The  problem  Is  then  to  find  the  function  w(x*y)  which  minimizes 
7 and  satisfies  the  prescribed  boundary  conditions  for  the  system  under 
consideration.  Direct  application  of  the  calculus  of  variations  to  min- 
imize 7 leads  to  the  analytical  problem;  that  la*  to  the  governing 
differential  equation  together  with  the  so-called  "natural'1  boundary  con- 
dltlona.  Instead  of  following  this  procedure,  Rltz  assumed  the  deflection* 
w(x,y),  ae  an  infinite  series  of  "admissible"  functions  and  adjusted  the 
coefficients  in  this  series  so  as  to  minimize  the  total  potential  energy* 
7.  To  be  admissible*  a function  must  satisfy  the  prescribed  boundary 
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conditions  of  the  problem.  Although  it  is  not  required , It  is  desirable 
to  choose  functions  which  also  satisfy  the  natural  boundary  conditions* 

In  order  to  be  assured  that  the  asstned  series  of  functions  can  accurately 
approximate  the  actual  deflection,  It  is  also  necessary  that  the  functions 
form  a complete  set.  If  now,  the  assuaed  deflection  function  is  sub- 
stituted into  the  integral  equation  for  the  total  potential  energy,  this 
energy  is  obtained  as  a quadratic  function  of  the  unknown  coefficients* 

The  variational  problem  is  thus  reduced  to  an  ordinary  minimum  problem* 

On  equating  to  zero  the  partial  derivative  of  V taken  with  respect  to 
each  coefficient,  there  is  obtained  an  infinite  system  of  linear  simul- 
taneous equations  which  define  the  values  of  these  coefficients.  An 
approximation  to  the  deflection  function  is  obtained  by  solving  a limited 
mmber  of  these  equations* 

5.  Potential  Energy  of  the  Loaded  Plate-Beam  System 

Consider  a plate  carrying  a uniformly  distributed  load  and 
supported  on  a rectangular  grid  of  beams  with  columns  at  their  inter- 
sections as  was  previously  described.  The  flexural  strain  energy  in  one 
panel  of  the  plate  is 

Uf*-  = JJA  {f  [w«+vv]  -»(■ -/*)[«<«  - w/y]]<JA  (D 

where  A is  the  region  within  the  panel  boundaries.  The  potential  energy 
of  the  uniform  load,  q,  on  the  plate  is 


W = - j JA  1 w JA 


(2) 
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and  the  flexural  strain  energy  of  a bean  of  length  L Is 


Consequently,  for  coordinate  axes  directed  along  the  beams,  the  total 
potential  energy  of  one  panel  of  the  loaded  plate  is 

V ” { f [w« + W^1 Z ~ D('  w****y  * wwJ)  d*dy 

♦fts  dx<,y  • 

The  deflection  is  zero  at  the  corner  of  each  panel  and,  frcn 
symmetry,  the  slope  of  the  deflection  surface  at  an  interior  panel 
boundary  in  the  direction  normal  to  the  boundary  is  also  zero*  Therefore, 
the  problem  reduces  to  that  of  finding  a symmetrical  function  v(x,y)  which 
minimizes  V,  given  by  equation  (4),  and  satisfies  the  imposed  boundary 
conditions, 

w (0,0)  = w(0,b)  = 0 

, v . , (5) 

wx(0,y)  - WyUfi)  - 0 

If  the  methods  of  the  calculus  of  variations  are  applied  to  equation  (4), 
one  finds  for  the  first  variation  of  the  potential  energy 

C Eq let  WXXXX (* P)  “ 2Pw yyy(*,0)]  Sw(XyO)  dx 

y)]Sw(0,y)  dy  (6) 

.L  *5 

+ \Q\o  [^74vv  " °|]  Sw(x,y)  dx  dy 


* 0 
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Since  8v(x,y)  is  an  arbitrary  function  which  need  satisfy  only  the  imposed 
boundary  conditions,  each  bracketed  quantity  must  vanish  identically. 
Rendering  V stationary  thus  implies  satisfaction  of  the  governing  differ- 
ential equation  and  of  the  natural  boundary  conditions^  On  applying  the 
variational  methods  to  the  second  term  in  the  integral  for  the  strain 
energy  of  the  plate,  it  is  found  that  the  boundary  conditions  alone  are 
sufficient  to  insure  that  the  variation  of  this  integral  vanish.  It 
thus  follows  that  this  tern  is  invariant  and  can  be  dropped  fraa  equation 
(4).  H.  Langhaar^  has  proven  that  this  ten  is  Invariant  when,  and 
only  when,  the  boundary  is  a geodesic  of  the  deformed  surface.  L.  E. 
Goodman  has  observed  ^ that  far  a uniformly  loaded  infinite  plate 
supported  on  either  (a)  a rectangular  gridwark  of  flexible  beams  support- 
ed at  their  inter  sections  by  column,  the  beams  in  either  direction  be- 
ing of  uniform  stiffness,  or,  (b)  a gridwark  of  flexible  beams  in  the 
shape  of  a regular  polygon  also  supported  at  the  corners  by  columns,  the 
beams  form  a boundary  along  which  the  deflection  of  the  deformed  surface 
ia  a geodesic. 
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III.  V.  J.  DUNCAN'S  S-FUNCTIONS 


6.  Definitions  and  Boundary  Properties  of  the  S-Functions 

A set  of  simple  polynomial  functions  with  properties  rery 

(7) 

suitable  to  the  present  problem  has  been  constructed  by  V.  J.  Duncan  • 
Since  only  the  symmetric  functions  have  been  used  in  the  present  work, 
they  will  be  the  only  ones  discussed*  The  general  symmetric  function 
can  be  defined  as 

c*  rT\  _ /4m  + 1 d2m  2 f ~zmr , . . 

Sm(5)  ('-»  J.  , (7) 


or, 

_ (4m-|)P,wtffet-0  -Zbm+OHJlZl -fl  +(4mt3)g„.2(2;-l) 

' 4(4m-lX4m+3)/4^TT  ’ (8) 

in  which  Pn  is  the  Legendre  polynomial  of  the  m-th  degree*  The  corres- 
ponding sequence  of  anti-syanetrlc  functions  can  be  obtained  from  the 
sequence  for  S (?)  by  replacing  m with  (m+l/2)  • The  polynomial  forms 
of  the  first  three  symmetric  functions  are 

S,(f)=  £r(sz-Zf3  + g*) 

S2(?)  = jr(3gz-20g3  + 45  5*-  4-2 f1-*- 14  f*) 

S3(g)  = ^-(gl-l4-g3  + 70f*-\6Bgs+ZI0i*-  I32S r+33fe). 

The  first  and  second  differential  coefficients  of  the  functions  are 


+ ~ ljw-,(Z»-l) 

2 j4m+\ 1 


(10) 
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and 


Sm(S)  = \/4m+l  1^(2? -l)  . ^ 

These  functions  were  derived  for  application  to  doubly  built-in  be  ana  v 
consequently  they  satisfy  the  following  boundary  conditions. 


5mC0)  = Sm(i)  = 0 

s ‘jo)  = s;(d  = o . 


(12) 


Numerical  values  of  the  first  four  functions  and  their 
second  derivatives  are  presented  in  Table  1.  nils  table  was  prepared 
with  the  aid  of  the  British  Association  Tables  of  Legendre  Polynomials 
(8). 


7.  Integral  Properties  of  the  S-Functions 

Lx  applying  sequences  of  functions  to  engineering  problems, 
it  is  necessary  to  evaluate  various  Integrals  which  contain  the  functions 
and  their  derivatives*  The  following  integrals  comprise  those  which 
occurred  in  solving  the  problem  under  discussion*  These  integrals  have 
been  evaluated  using  equations  (8)  and  (11)  and  well  known  integral 
properties  of  the  Legendre  polynomial  functions  (9)* 


f 

’o 


j: 


SJl)  df 

= 0 

- JT 
'Jo 

m 1 

m * 1 

(13) 

s;(«d? 

= 0 

all  m 

(U) 

S;«  S„*(f)  dS  = 0 

n 

(15) 

= 1 

n - m 

13. 


f Sm(5)S»df  = 0 


I 


4(4  m +3)  /(4  m + 1 ) (4  m +5) 

-I 

Z (4  m - l)(4m  +3) 

1 

4(4  m-l)/ (4m  +l)(4m~  3) 


n & m , n m ± I 
n = m + l 
n = m 
n -m  - I 


(16) 


r sjf)^)  df  = 


n*m  , n#m  ± I , n -m±2 


8 (4  m - 1 )(4  m -3X4  m + 3)  ( 4 m 5) 

-I 


4(4  m-  l)(4m+3Xlm"5)/(4m-fl)(4m-3/ 


I 


^(4m-lXim-3)(4m-5)/(4m+|)(4m-7) 


n=  m 
n - m- 1 
n = m— 2 


(17) 


jXod*  -f 

= 8(4m  -l)(4m-3;(4m+3)(4m+5;/4W7f 

- 3(4m-3)(4m  *3)gm>l(2?-lJ  +3  (4  m -l)(4m  +5)Ptm-,(2S-l) 

- (4m+3X4m+s)rim.3teT-0j  t m >1 


2 /4m  + I 


(19) 


Numerical  values  for  sane  definite  integrals  of  the  fonts  of 


equations  (16),  (17),  (13)  and  (19)  are  given  in  Table  2, 
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17.  SOLUTION  QF  THE  PROBIZM  USUG  THE  S-FUBCTICKS 


8.  The  Equations  for  the  Deflection  Function  Coefficients 

It  is  dosirod  to  asstme  the  deflection  function  w(x,y)  as 
a soquanoo  of  symmetrical  functions  which  ara  capable  of  approximating 
tha  daflactod  shape  of  ths  plat*  as  closely  as  may  be  desired.  Each 
function  of  this  sequence  must  satisfy  the  prescribed  boundary  conditions 
given  by  equation  (5 ) • These  conditions  are  fulfilled  by  asststing  a 
deflection  function  in  the  following  fora. 


^qVd  = 2j  5m(f)  SJf)  + $n(l[)  . (20) 

The  coefficients*  <Xrm,  and  f are  to  be  determined  such  that  the 
total  potential  energy  7 will  be  a minimum.  Omitting  the  second  term 
from  the  first  integral  of  equation  (4)  and  writing  the  equation  in 
dimensionless  farm*  the  expression  for  7 becomes. 


(21) 


v = °*c  1 £ {2?K}+'WcJ]2'  W df 

The  energy  7 is  obtained  as  a quadratic  function  of  the  coefficients 

and  Xn  by  substituting  equation  (20)  into  equation  (21)*  that  is* 

* 2 i I {|  S*4sw®5,,(7)  2fa)J 

+%$mS"jg)  + yBS„*(7))2df  dr, 

~ Jo  Jo  HP  07  (22) 

+ ££{$  t + d7>  ■ 
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Let  r be  any  value  of  a and  a be  any  value  of  n.  Oh  the  assumption  that 
the  order  of  suoaatlons  integrations  can  be  interchanged  and  using 
the  properties  of  the  S-functlons  given  by  equations  (13) , (14)  and  (15 )» 
then. 


The  coefficients  are  to  be  chosen  such  that  V will  be  a 
nlnlwm;  thus  for  the  coefficients  and 


il  = iZ  = o «) 

a** 


which  load  to  the  relationships. 


„ _ yy  (i -<*..) 

" 60  (I  + \a/c) 

a _ yy  ^(n  i 

"m  60  (l  + Ve) 


m * I 


(25) 


y _ yy  (c^-ofr) 

J|  60  (l  + cXb) 

y _ yF  o<M 

*"  60  (i  + c>k) 


n*l 
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Ckx  substituting  these  values  of  ft.®*4  •«  into  equation  (23)*  the 
potential  energy  Is  obtained  as  a function  of  the  coefficients  ofmn 
alone.  Making  this  substitution  and  dropping  a constant  tent,  coo  obtains* 


_ V 

+ 


= ????  ?^£s;srdf  f S„Ss"d»} 


s^f's.s.dT,  45? 


(26) 


» f 

\440{\+K£)  2 


c< 


mi 


1440c4(|+cXb) 


§«S 


o<m  f ! + _ , — i — -1  + °*il (! — 

1440  LO +*«/£)  c*0+cAj j 720  (l+A^yO+c?^) 


Using  this  equation  and  the  ntsterlcal  ralues  of  the  Integrals  given  In 
Table  2,  the  matrix  of  Table  3 far  the  total  potential  energy  has  been 
obtained.  The  total  energy  Is  equal  to  the  Uniting  sun  of  the  entries 
In  the  table  multiplied  by  the  coefficients  at  the  heads  of  their  res- 
pect ire  rows  and  colmna,  It  should  be  noted  that  the  natrlx  Is 
syetrlcal  about  Its  principal  diagonal. 

The  energy  7 is  minimized  by  equating  the  partial  derivative 
with  respect  to  each  coefficient  o<Mto  aero.  Consequently*  the  linear 
simultaneous  equations  defining  these  coefficients  are  obtained  by 
equating  each  row  of  Table  3 to  aero  and  cancelling  the  factor  n2"  at 
the  top  of  the  column  of  constants.  The  properties  of  the  S-functlona 
given  by  equations  (13)*  (14)  and  (15)  have  resulted  in  the  staple  re- 
lationships of  equations  (25)  between  the  coefficients  ^ and  Xj,  and  the 
coefficients  tfmn*  and  thereby  have  effected  a considerable  reduction  In 
the  number  of  simultaneous  equations  to  be  solved. 
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9.  Moments  and  Afar  age  Momenta  in  the  Plata  and  Moownta  In  the  Beams* 
Momenta  In  the  plat#  and  beans  are  found  tram  the  deflection 
function  v In  the  uaual  manner*  The  plate  moments  are 


in  which  the  ■aero"  superscripts  are  uaed  to  indicate  that  the  moment  is 
given  for  a value  of  Poisson1 a ratio  equal  to  aero.  Equations  (27)  are 
very  useful  for  computing  Momenta  at  a point  for  any  desirable  value  of 
Poisson's  ratio*  If  Poisson's  ratio  ia  taken  aa  zero,  the  aonenta  are 

= <*mn  S*f5)$,fy)  - <jf a*  j*  ^ Si  Cl) 

(28) 

My  =-  ^ S.  °<mnSm(f)  5n(7f)  - ^ }fn  Sn(? }) 

In  addition  to  a knowledge  of  aonenta  at  particular  points  in 
the  plate,  it  ia  often  desirable  to  know  the  naxlatn  average  moments  in 
what  are  commonly  called  the  ■middle " and  "column”  stripe  of  the  panel. 
The  term  "column  atrip"  la  here  used  to  designate  the  portion  of  the 
panel  over  a beam  and  its  supporting  column  and  having  a width  equal  to 
half  of  the  panel  width;  the  term  "middle  strip”  is  used  to  designate 
the  middle  portion  of  the  panel,  one-half  the  width  of  the  panel.  The 
average  moments  on  the  center  and  end  sections  of  the  various  ooltam 
and  middle  stripe  are  defined  as  follows t 
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M,  - M®  +/N,  « 4 J*  ^(£-,7)  cty 

M*=M£+/^  = 4 + 4/  J^MyO^) <hj 

• ^ 

M3  = Ml  +yu  N3  = 2 JfKftfl)  d«f  + 2^  M i&T,)  J, 
m4  - m;  v N4  = 2 J^Mx*(i,if)  ^ + 2/  d>? 

ms  - m;  +/*ns  - 4jf m;(i  i)  d?  + v JoV(?.i)  d ? 

Ms  - m;  +*k,  ■ )df  + v (*m;M  j? 
Mr  = M’7  v*N7  * 2J^Myfe^dr  + 2/ 

M8  = Me  +*  Ne  = 2j*Mp(?,i)dS  +2/J*M;(f,l)df 


(28) 


The  mzaerical  Talus 8 of  the  integrals  encotmtsrsd  aftsr  substituting 
equations  (27)  into  aquations  (28)  are  given  in  Table  2.  Since  the 
deflection  surface  is  sywaetrical  and  has  zero  slope  at  the  panel  bound- 
aries, the  tens  (H^+H^),  (H^I^),  (K^+fc^)  and  (I^+Hg ) must  each  be 
zero.  The  average  momenta  on  a full  central  or  edge  section  are  there- 
fore independent  of  Poisson's  ratio.  The  locations  of  the  lines  along 
which  the  moments  are  averaged  are  shown  in  Figure  1. 
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FIG.  1 LOCATION  OF  AVERAGE  MOMENTS 

The  moments  in  the  beau  can  be  determined  from  the  coefficients 
and  or  they  can  be  obtained  directly  from  the  plate  moments* 

Mq  = ~1*\,  Jft.  ° 

(29) 

Mb  = - 2 t„  S^)  = a c Ak  MJM  . 

10.  Reduction  to  an  Exact  Solution  when  = 1* 

It  is  evident  from  Table  3 that  the  coefficients  o^vill  all 

be  sero  If 

A„  Ak  - l , (30) 

consequently.  It  follows  from  equations  (25)  that  the  only  non-zero 


coefficients  are  §,  and  • On  substituting  ths  values  of  tbs  co- 
efficients into  aquation  (20),  ths  ds flection  is  found  to  bo 
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w = & ^£1  r 5,(f)  + °4  Si(7l)  1 

60  D l l + A«/c  l+cAbJ  > 

or,  replacing  S-^  by  its  polynomial  expansion  given  in  equations  (9)» 


r*0  -Si1 


24 D L l+>vfc 


I + c Ab 


] . 


and. 


For  Poisson’s  ratio  equal  to  aero,  the  moments  in  the  plate 


" ,AV*)  (^7  . 

It  is  interesting  that,  far  this  case,  the  moment  is  equal  to  the 

moment  in  a clamped-ended  beam  of  unit  width  and  of  length  equal  to  a, 

under  a uniformly  distributed  load  of  per  unit  length}  Kf 

(l  + A«/c;  * 

is  equal  to  the  moment  in  a similar  beam  but  of  length  equal  to  b,  and 

<? 

tinder  a load  of  (|  + c X J P®r  t»it  length#  Therefore,  plates 
supported  on  beams,  with  relative  rigidities  such  that  W~  1* 
function  much  like  independent  prismatic  beam  elements}  the  load  is 
distributed  in  the  ratios  ls(lec  \)  and  1:  (l+A#/c)  to  the  short  and 
long  elements  respectively#  There  are  no  twisting  moments  in  such 
plates# 


The  shearing  forces  in  the  plate  are 


The  supporting  beans  parallel  to  the  x-axis  (long  beams)  are  therefore 
subjected  to  a uniform  load  of  per  unit  of  length  while  the  load 

I +A«/c 

QQC 

on  those  parallel  to  the  y-axia  is  ~ 


per  unit  of  length. 
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V.  TABI£S  AND  GRAPHS  OF  NUMERICAL  SOLUTIONS 
11*  Tables  of  Results 

Numerical  solutions  hare  been  made  for  about  forty  eases , 
covering  three  different  panel  side-ratios,  c — 1*0,  c = 0.8,  and 
c = 0.5,  and  embracing  a vide  range  of  beam-plate  rigidities,  Aa  and  A^ 

In  the  ease  of  square  panels,  the  rigidities,  A#  and  A^,  were  chosen  so 
as  to  represent  all  combinations  of  values  from  zero  to  Infinity.  For 
the  other  panel  side-ratios,  it  vas  assumed  either  that  all  beams  were 
of  equal  rigidity  (i.e.  Ab  = A^)  or  that  all  beans  had  equal  moments  of 
inertia  (i.e.  ^ = I#  or  Afe  = \/c).  A tabulated  stnmary  of  the 
solutions  is  presented  in  Table  4. 

La  general,  the  results  presented  in  Table  4 are  those  obtained 
Turing  nine  of  the  coefficients  o<m;  that  is,  m and  n were  taken  over  the 
range  from  1 to  3.  Each  case  thus  involved  the  solution  of  a set  of  nine 
simultaneous  equations.  In  addition,  each  case  vas  solved  using  only 
m = n = 1 and  again  using  m = 1,2  and  n = 1,2.  A sequence  of  approximate 
solutions  vas  thus  obtained;  the  rate  of  convergence  of  this  sequence 
could  be  estimated  by  inspection.  The  accuracy  of  the  results  presented 
in  Table  4 ia  indicated  by  the  underlining  of  digits  vhich  differ  ty 
more  than  five  units  in  comparison  with  a solution  using  fever  terms  in 
the  series.  For  example,  solutions  obtained  using  nine  coefficients  o^ 
vere  compared  to  the  solutions  obtained  using  only  four  such  coefficients. 
The  solutions  for  square  symmetrical  panels  vere  obtained  using  m = 1,2, 
3,4  and  n = 1,2, 3,4;  these  oases  required  the  solution  of  sets  of  ten 
simultaneous  equations.  The  results  for  the  limiting  case  of  a plate 
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supported  on  coltn&s  without  connecting  beams  were  obtained  from  the 
exact  Infinite  series  solution  which  is  siamarlzed  in  the  Appendix* 

12*  Graphs  of  Results 

The  results  are  presented  in  graphical  fora  in  Figures  2 to 
42*  In  order  to  plot  the  relative  bean  rigidity,  A#»  over  its  full 
range  tram  aero  to  infinity,  a modified  abscissa  scale  has  been  used* 

The  abscissa  scales  are  divided  in  the  ratio  (2+A^),  therefore 
^ = 2*0  falls  at  the  center  of  the  scale,  = 1*0  falls  at  the  1/3 
point  and  so  on*  Figures  2 to  15  give  deflections,  moments  and  average 
moments  far  a square  panel  far  all  combinations  of  be  am- plate  rigidities* 
Figures  16  to  42  give  these  deflections  and  moments  far  panels  with  side 
ratios  of  c = 1.0,  0.8  and  0*5,  for  cases  in  which  X^  = Xa  or  1^  = I#. 
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VI.  ACCURACY  Of  THE  SOLUTIGHS 

13*  Convergence  of  the  Sequence  of  Approximate  Solutions 

Cbe  of  the  most  undesirable  features  of  the  Ritz  procedure 
is  that  it  does  not  include  a method  or  general  principal  for  determin- 
ing rigorously  the  accuracy  of  the  results  (10,11)*  However,  some  in- 
dication of  the  accuracy  can  be  obtained  by  inspecting  the  convergence 
of  solutions  in  which  increasingly  larger  mttbers  of  terms  of  the  approx- 
imating series  were  used*  As  has  been  previously  described,  a sequence 
of  three  such  solutions  was  generally  obtained  for  each  case  considered 
in  the  present  work*  Because  of  the  simplifications  Implied  by  symmetry, 
it  was  feasible  to  obtain  a sequence  of  four  approximate  solutions  for 
cases  of  square  panels  in  which  These  solutions  are  simarized 

in  Table  5 end  in  the  first  page  of  Table  6* 

The  convergence  of  the  sequence  of  solutions  for  the  clamped 
plate  - oo)  is  very  good,  however,  the  convergence  becomes  even 

better  as  ^ and  Afe  approach  unity.  In  fact,  as  has  been  previously 
discussed,  the  solution  becomes  exact  when  A^A^  = 1-  As  Aa  and  A^  are 
decreased  further,  the  convergence  of  the  solutions  becomes  poorer*  The 
solutions  show  the  poorest  convergence  when  Aa  and  are  both  zero; 
the  results  of  these  eases  are  presented  in  Table  6*  From  this  table, 
it  is  apparent  that  the  convergence  is  better  when  the  panels  are  square 
than  when  they  are  long  and  narrow* 
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lA.  Limiting  Cases ; Plate  Supported  on  Columns  Only  and  the  Clamped 
Plate 

The  bending  of  a uniformly  loaded  plat*?  supported  by  rows  of 

(3  1 2) 

equidistant  columns  has  been  discussed  by  several  authors  v ’ • The 

analytical  solution  of  this  problem  leads  to  the  deflection  function  ex- 
pressed as  an  infinite  aeries  of  trigonometric  and  hyperbolic  functions. 
This  formula  for  the  deflection,  v,  as  well  as  formulae  for  moments  and 
various  average  moments  are  summarized  in  the  Appendix*  Numerical  re- 
sults obtained  using  these  formulae  are  summarized  in  Table  6 in  the 
columns  headed  "Exact  Solution" • 

The  agreement  between  the  energy  solutions  and  the  exact 
solutions  is,  in  general,  fairly  good*  Because  of  the  assumption  that 
the  columns  provide  only  concentrated  point  reactions,  the  exact  sol- 
utions shows  infinite  plate  moments  at  the  panel  corner" - Therefore,  the 
moments  obtained  using  the  energy  method  with  a finite  number  of  tezms 
cannot  agree  with  those  of  the  exact  solution  at  these  points*  Aside 
from  this,  the  least  accurate  moments  obtained  using  the  energy  method 
are  the  negative  momenta  at  the  middle  of  the  long  sides  of  the  panels; 
these  moments,  however,  are  generally  fairly  small* 

Another  Important  solution  obtained  as  a limiting  case  in  the 
present  work  is  the  solution  to  the  problem  of  a uniformly  loaded  clamped 
plate.  An  exact  analytic  solution  has  not  been  obtained  for  this  problem 
so  that  exact  results  are  not  available  for  comparison*  However,  certain 
deflections  and  moments  have  been  computed  by  other  investigators  using 
various  approximate  methods  (see  for  instance  references  (3)  and  (13)  to 
(21)).  Usually  the  method  of  Ritz  (described  herein)  or  the  methods 
developed  by  Hencky  ^)  or  Timoshenko  ^ are  used. 
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Son*  results  obtained  by  investigators  using  Hencky's 
analysis  are  presented  in  Table  7;  these  are  virtually  in  agreement 
with  results  from  the  present  investigation  which  are  also  presented 
in  this  table.  In  Hencky's  analysis,  the  deflection  function  consists 
of  a parabolic  tern  plus  a correction  function  in  the  fora  of  a series 
of  trigonometric  and  hyperbolic  functions.  This  deflection  function 
satisfies  the  differential  equation  (Dv^w  = q)  exactly;  the  co- 
efficients in  the  series  are  adjusted  so  that  the  boundary  conditions 
are  approximately  satisfied.  The  numerical  values  given  in  Table  7 
for  Hencky's  analysis  were  computed  by  Siess  and  Newmark  from  co- 
efficients tabulated  by  Vojtaszak  for  the  terms  n = 1,  3*  5 ....27 
of  the  series.  Average  moments  on  an  edge  were  computed  by  integration, 
however,  average  moments  on  Interior  sections  were  obtained  by  the 
application  of  Simpson » s one-third  rule  to  moments  at  the  one-tenth 
points. 

15 . Finite  Difference  Solution* 

The  calculus  of  finite  differences  offers  sn  alternate  approx- 
imate numerical  procedure  for  solving  the  general  problem  of  a plate 
continuous  over  flexible  beams.  A number  of  such  solutions  have  been 
made  at  the  University  of  Illinois  by  various  persons  working  under  the 
direction  of  Professors  N.  M.  Newmark  and  C.  P.  Siess.  These  solutions 
were  made  using  a difference  net  which  divided  each  span  into  eight  equal 
parts.  Each  panel  vas  therefore  divided  into  64  rectangular  elements. 
With  this  difference  net  and  noting  the  simplifications  resulting  from 
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synnetry,  It  is  found  that  each  case  requires  the  solution  of  24 
simultaneous  linear  equations.  For  square  plates  supported  on  beans 
which  are  all  of  the  sane  rigidity,  the  ntnber  of  equations  reduces  to 
14. 

Some  of  the  deflections  and  sonants  obtained  by  these 
difference  equation  solutions  are  presented  in  Tables  7 and  8.  It  is 
Interesting  to  note  that  the  "S-function”  solutions  using  only  four  co- 
efficients (four  simultaneous  equations)  are  generally  rare 
accurate  than  the  difference  solutions  using  24  slnultaneous  equations*. 
It  is  also  interesting  that  the  difference  equations  determine  the 
noaents  rare  accurately  than  they  do  the  deflections. 
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VII.  CONCLUSION 


16.  Sumaary  and  Concluding  Bsmarks 

Tha  present  investigation  has  resulted  in  a relatively 
s in pie  approximate  analytical  solution  to  the  problem  of  an  interior 
panel  of  a uniformly  loaded  plate  which  is  continuous  over  a rect- 
angular grid  of  flexible  beams  supported  at  their  intersections  by 
columns.  The  relative  simplicity  of  the  analysis  made  it  possible  to 
present  fairly  accurate  and  complete  results  for  panels  with  various 
side  ratios  and  various  ratios  of  beam  rigidities  to  plate  rigidities. 
These  results  are  given  in  Table  4 and  Figures  2 to  42. 

Although  the  results  presented  are  those  for  a plate  uniform- 
ly loaded  over  its  full  area,  they  can  be  used  to  obtain  solutions  for 
several  other  cases  of  loading.  In  particular,  solutions  can  be 
obtained  for  the  case  in  which  only  alternate  panels  are  uniformly  load- 
ed and  for  the  case  in  which  only  alternate  rows  of  panels  are  uniformly 
loaded.  The  alternate  panel  loading  can  be  obtained  by  superposing  a 
uniform  load  q/2  with  another  uniform  load  q/2  which  is  alternately 
positive  and  negative,  changing  sign  at  each  panel  boundary.  TJhder  this 
latter  loading,  and  assuming  that  the  torsional  restraint  of  the  beams 
can  be  neglected,  each  panel  is  in  the  same  condition  as  a simply 
supported  plate  of  panel  size.  Therefore,  when  alternate  panels  are  load- 
ed, the  solution  can  be  obtained  by  superposing  the  present  solution  with 
that  for  a simply  supported  plate.  Similarly,  the  uniform  loading  of 
alternate  rows  of  panels  is  equivalent  to  superposing  a uniform  load  q/2 
with  a uniform  load  q/2  which  changes  sign  from  row  to  row  of  panels  but 
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which  has  the  same  sign  in  any  one  row*  Under  this  latter  loading  each 
row  of  panels  is  essentially  a simply  supported  strip.  Otoe  representative 
panel  from  this  strip  can  he  solved  by  Levy’s  method  ind  oan  be 

superposed  with  the  solution  given  herein  to  obtain  the  final  solution. 

One  of  the  most  important  types  of  structures  similar  to  the 
plate-beam  system  described  is  the  type  of  reinforced  concrete  floor 
known  as  the  two-way  slab.  The  foregoing  solution  is  based  on  a ntnber 
of  asstnptions  which  are  not  completely  satisfied  by  this  type  of  struc- 
ture. The  results  presented  therefore  require  some  modification  before 
they  can  be  applied  in  the  design  or  analysis  of  such  two-way  slabs. 

The  purpose  of  this  investigation  was  not  the  actual  development  of  an 
analysis  of  reinforced  concrete  slabs.  However , analytical  solutions, 
such  as  those  presented,  are  a basic  step  toward  the  ultimate  objective 
of  a rational  design  procedure. 
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Sunn  ary  of  the  Szact  Solution  for  tha  Casa  \ = \ = 0 


The  following  summary  of  formulae  far  deflections,  momenta, 
and  average  sonants  in  a uniformly  loadad  plate  supported  by  rows  of 
equidistant  coltmns  is  based  on  tha  analytical  solution  presented  in 

(o\ 

section  46  of  S.  Timoshenko  "Theory  of  Plates  and  Shells"  ' • The 
notation  used  follows  that  of  the  previous  part  of  this  work,  however, 
the  origin  of  coordinates  is  taken  at  the  center  of  a panel.  The  de- 
flection of  the  panel  of  the  plate  is 
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where 


in  alternate  form  of  the  deflection  function  is 


obtained  by  interchanging  with  y^  and  a with  b. 
The  bending  moments  are* 
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The  ererage  moments  in  the  coltnm  and  middle  stripe  ares 
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and  the  noment  coefficients  for  correcting  these  average  moments  for 
Boissons  ratio  are: 


= Jk-=  _s!  + 

<jq*  qqI 2-  32  8 


Jk  - Jk  - + _s! 

“ qq*  "32  8 


(-l)fr  tanh  (**Jz 
. ©<m  cosh  cX^/2 


fonh  oWg 
©ha  cosh  cK*/i 


Ny  _ _ | ^ 1 (-ift  tflnh  < Xm/Zc* 

q?  32  8 ££  o<m/c*>  eosh<*w/2e* 


- JSkL- Jk  = 
qo*  qo* 


I , _L  V fonh  o<m Ac* 

"3?  8 <*m/c*  cosh  <**x/2.cx 


TABUS  1.  THE  FIRST  FOUR  S - FUHCTIQBS  AHD 
THEIR  SECCHD  DERIVATIVES 


s 

^(f) 

0.00 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

♦0.000000 

♦0.002523 

♦0.009056 

♦0.018175 

♦0.028622 

♦0.Q39306 

♦O.Q49305 

♦0.057865 

♦0.064399 

♦0.068487 

♦0.069877 

♦2.236068 

♦1.598789 

♦1.028591 

♦0.525476 

♦0.089443 

-0.279508 

-0.581378 

-0.816165 

-0.983870 

-1.084493 

-1.118CB4 


♦0.000000 
♦0.002634 
♦0.007047 
♦0.009676 
♦0.009728 
♦0.006592 
♦0.001323 
-0.004787 
-O.0QL0368 
-o.  014242 
-0.015625 


s;cf) 


♦3.000000 

♦0.623813 

-0.699000 

-1.236187 

-1.224000 

-0.867187 

-0.339000 

♦0.218813 

♦0.696000 

♦1.013813 

♦1.125000 


•00 
.95 
.90 
.85 
•80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 


? 

S3  (8) 

0.00 

♦0.000000 

0.05 

♦0.002052 

0.10 

♦0.002735 

0.15 

♦0.000189 

0.20 

-0.003471 

0.25 

-0.005694 

0.30 

-0.0051U 

0.35 

-0.002057 

0.40 

♦0.002160 

0*45 

+0.005683 

0.50 

♦0.007042 

♦3.605551 
-0.869530 
-1 .412641 
-0.451726 
♦0.620501 
♦1.165466 
♦1.055U4 
♦0.465769 
-0.290521 
-0.897167 
-1.126735 


s4(f) 

».  . 

S4(f) 

♦0.000000 

♦0.001149 

-0.000504 

-0.002451 

-0.001618 

♦0.001239 

♦O.OQ3423 

♦0.0Q3054 

♦0.000372 

-0.002700 

-0.004026 

♦4.123106 

-1.689178 

-0.068672 

♦1.264574 

♦0.875497 

-0.303621 

-1.100666 

-0.985730 

-0.163130 

♦0.743481 

♦1.127412 
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TABI£  2.  NUMERICAL  VALDES  OF  VARI0D5  INTEGRALS  INVOLVING 
THE  S - FUND TICKS  AND  THEIR  DERIVATIVES 


(a)  Valuta  of  (10H  S_(!)S_(f)df 


♦1964*1270 


2 

3 

-161.3325 

♦74.92508 

♦11.18639 

-20.01083 

♦14*85443 

♦2.355630 

-5*364217 

♦4*819742 


Cb)  Valuta  of  (10)3J  S*(f)Sn(f)df 


-23.80952 
♦ 5.323971 
0 
0 


♦5.323971 

-6.493506 

♦2.101137 

0 


♦2.101137 


♦1.121121 

-1.754386 


(o)  Mtan  valuta  of  8^(1)  and  S (J)  par  quartar  span 


♦0.01543120 

♦0.00659180 

-0.00024757 

-0.00060161 


2L  SJW 


♦0.05910440 

-0.00659180 

♦0.00024757 

♦0.00060161 


♦0.8385255 

-0.3515625 

-0.C739420 

♦0.2381909 


TABUS  3.  MATRIX  PGR  THE  TOTAL 


TABIE  4.  EKFIZCTICBS  AND  MOMENTS  15  TEE  PLATE  AID  BEAKS  FOR 
VARIOUS  BEAM  STIFFNESSES  AND  VARIOUS  PANEL  SID3  RATIOS 
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..0 

Units 

06 

03 

exS 

cx3 

c*d 

5.0 

2.0 

1.0 

§ 

exit 

ex23 

Sh 

Table  5 

♦0.2453540 

♦0.0502558 

♦0.0452006 

+0.0347526 

♦o. 0014119 

♦0.0009518 

♦0.0006044 

♦0.0003340 

-0.0062936 

♦0.2069246 

♦0.0484327 

*0.0386811 

♦0.0369009 

♦0.0020701 

♦0.0010036 

♦0.0040019 

♦0.0Q31922 

-0.0025925 

♦0.1642342 

♦0.0448687 

♦0.0313020 

♦0.0374465 

♦0.0028536 

♦0.0010158 

♦0.0074132 

♦0.0057173 

+0.0011657 

V (£,£) 

w(0,J) 

w(i*0) 

0.001265 

0 

0 

0.001435 

0.0003324 

0 

0.001622 

0.0006977 

0 

0.001829 

0.001101 

0 

q.Vn 

qa‘/D 

q«VD 

Sit) 

-«(o,o) 

-19(0,0) 

0.01762 

0.01762 

0 

0 

0.01537 

0.02068 

0 

0.009512 

0.01292 
0; 02405 
0 

0.02033 

0.01020 

0.02777 

0 

0.  Q3212 

q-2 

3 

q* 

-My  (i*£) 
My(0#£) 

0 

0.05132 

0.05132 

0 

0 

0.05536 

O.Q4452 

0.005582 

0 

0.06481 

0.02906 

0.01829 

q*2 

*2 

q*2 

qa 

-5 

•5 

-MJ 

H* 


111 


5a  $ 

-jrfoj 


0.00410^ 

O.OQ3581 

0.003025 

0.002404 

0.0132S 

0.01205 

0.01050 

0.008702 

0.01505 

0.01315 

o.onoS 

0.008745 

0.04492 

0.Q3920 

0.Q3290 

0.02594 

O.OQ4107 

0.00904s 

0.01447 

0.02041 

0.01328 

0.0217S 

0.03102 

O.Q4140 

1 

0.01845 

0.02634 

0.04492 

O.Q4975 

0.055 n 

0.06lCg 

0.01467 

0.01678 

0.01910 

0.02166 

0 

O.OQ4Q41 

0.008475 

0.01336 

0.01467 

lit 

0.01074 

0.008471 

0 

0 

0 

0 

0.  Q3209 

0.03330 

0.03462 

0.Q3609 

0.05423 

0.05m 

O.O6I65 

fill  1 

0.03209 

0.02791 

0.02333 

0.04756 

0.03212 

qa2 

qa2 

£ 

q*2 

q*2 

£ 

q*2 

^2 

S2 

q*3 

q*3 

q*3 

qa 


TABLE  4*  ( continued) 


*12 

°122 

S31 

S23 

32 

o<33 


w(i,o) 


i!S 


♦0.1164395 

♦0.0382533 

♦0.0228Q39 

♦O.Q348641 

♦0.0036594 

♦0.0009488 

♦0.0103083 

♦0.0073828 

♦O.OQ45744 


0.002059 

0.001549 

0 

0.0G71T2 

0.03191 

0 

0.04594 

0 

0.0703S 

0.02022 

0.02551 

O.OOOL7 


0.002604 

0.002604 

0 

0 

0.04167 

0 

0.08333 

0 

0.08333 

0 

0.04167 


0.04167 

0.08333 

0.04167 

0.08333 

0.Q3125 

0.Q3125 

0 

0 

0.04167 

0.08333 

0 

0 


0.001622 

O.OOQ3484 

0.0003484 

0.01849 

0.01849 

0.01019 

0.01019 

0.005797 

0.04857 

0.04857 

0.005797 

0.008732 

0.02072 

0.01663 

0.04400 

0.008732 

0.02C72 

0.01663 

0.G4400 

0.01492 

0.004226 

0.01492 

0.004226 

0.02899 


♦0.1640118 

♦0.0388454 

♦0.Q344938 

♦0.Q341673 

♦0.0018214 

♦0.0012804 

♦0.0051238 

♦0.0046578 

-O.OOQ3345 


0.001827 

0.00C7331 

0.0003661 

0.01605 

0.02191 

0.01094 

0.02162 

0.006032 

0.05294 

0.041(g 

0.01213 

0.008352 

0.01962 

0.01460 

0.03762 

0.01436 
0.Q305 


O.Q493& 

0.01726 

0.006880 

0.01287 

0.004430 

O.03Q19 

0.05470 

0.Q242S 

0.04332 


TABXE  4*  (continued) 


<*£ 

2a 

S22 

£31 

o<23 

S32 

33 


&8 


o = o 


+0.1161059 

♦0.0321086 

♦0.0249051 

♦0.0301154 

♦0.0021838 

♦0.0010926 

+0.0067431 

♦0.0056569 

♦0.0021525 


0.002058 

0.001160 

0.0003860 

0.01335 

0.02572 

0.01179 

Oi 03482 

0.0063CJ2 
0.05784 


♦0.0619847 

♦0.02C6235 

♦0.0136689 

+0.0206320 

♦0.0020724 

♦0.0007295 

♦0.0063941 

♦0.0048797 

♦0.0Q32376 


0.002315 

0.001636 

0.0004084 

0.01Q34 

0.02996 

0.01276 

0.05032 

0.006604 

0.06331 

0.0237S 

0.02657 

0.0C7482 

0.01632 

0.009802 

0.02271 

0.02732 
0.0543 


r.w.vi.T?i 


0.06192 

0.02277 

0.01973 

0.006052 

0.004916 

0.Q3302 

0.06372 

0.01329 

0.02520 


-0.C727797 

-0.0418028 

-0.0130433 

-0.0487697 

-0.0116478 

-0.0017080 

-0.0304619 

-0.0166756 

-0.0200419 


0.002932 

0.002772 

0.0004639 

0.003921 

0.04007 


0.007336 

0.07631 

0.0Q3<g0 

0,04344 

0.006224 

0.009902 

0.0Q3682 

0.0C4Cgg 

0.04271 

0.06861 

0.04062 

0.07892 

0,02974 

O.O3294 

0.002012 

0.005547 

0.03668 

0.07632 

0 

0 


0.002058 

0.0007726 

0.0007726 

O.Q1954 

0.01954 

0.02345 

0.02345 

0.01265 

0.04534 

0.04534 

0.01265 

0.01422 

0.02992 

0.01852 

0.04293 

0.01422 

0.02992 

0.01852 

0.04293 

0.01523 

0.009332 

0.01523 

0.009332 

0.02529 

0.0(691 

0.08629 

0.0(691 


q4/B 

q«/D 


I 


T1EZ£  4*  (continued) 


Unite 


♦0.0618319 

♦0.0173684 

♦0.0151713 

♦0.0185623 

♦0.0012795 

♦0.0008993 

♦0.0046452 

♦0.0042376 

♦0.0021341 


0.002315 

0.001226 

0.0008173 

0.01688 

0.02342 

0.02546 

0.03802 

0.01323 

0.05015 

0.03698 

0.01988 

0.01407 

0.02906 

0.01634 

0.03579 

0.02068 
0.0417 


0.04886 

0.01787 

0.01477 

0.01296 

0.009840 


0.002604 

0.001736 

0.0008681 

0.01389 

0.02778 

0.02778 

0.05556 

0.01389 

0.05556 

0.02778 

0.02778 

0.01389 

0.02778 

0.01389 

0.02778 

0.02778 

0.05556 

0.02778 

0.05556 

0.02083 

0.02083 

0.01042 

0.01042 


-0.1578685 

-0.0924806 

-O.Q449433 

-0.1149632 

-0.0266033 

-0.0050624 

-0.0715104 

-0.0412662 

-O.Q483329 


0.003311 

0.002968 

0.0009968 


0.002604 

0.001302 

0.001302 

0.02083 

0.02083 

0.04167 

0.04167 

0.02083 

0.04167 

0.04167 

0.02083 

0.02063 

0.04167 

0.02083 

0.04167 

0.02063 

0.04167 

0.02063 

0.04167 

0.01562 

0.01562 

0.01562 

0.01562 


U.U0O4O 

0.05091 

0.01983 

0.Q38C2 


U.U*f70 

0.05556 

0.01389 

0.02778 


u,u<:ue^ 

0.04167 

0.02083 

0.04167 


42. 


ZABLE  4»  (continued) 


w(i,|) 

wfe.fi 

w(i,0) 


-HyClfO) 

(Otl/ 

M? (0,i) 


iiS 


-0.0714177 

-0.0246824 

-0.0210902 

-0.0311203 

-0.0028365 

-0.0019368 

-0.0106518 

-0.0095379 

-0.0069255 


0.002933 

0.001851 

0.001389 

0.01790 

0.02530 

0.04595 

0.0616^ 

0.02192 


0.04121 

0.01850 

0.03355 

0.02829 

0.05662 

0.02589 

0.04836 

0.01865 

0.02210 

0.01310 

0.01660 

0.02192 

0.04595 

0.01457 

0.03081 


-0.2592456 

-0.1559400. 

-0.0867370 

-0.2095022 

-0.0467405 

-0.0119897 

-0.1298752 

-0.0803875 

-0.0910728 


0.003754 

0.003201 

0.001612 

0.01068 

0.03616 

0.05*12 

0.l2g£ 

0.02462 

0.0596^ 

0.01042 

O.Q4722 

0.02130 

0.03625 

0.01270 

0.01402 

0.04522 

0.1015 

0.03806 

0.065M 

0.02607 

0.0371g 

0.0069% 

0.01909 


0.003312 

0.001985 

0.001985 

0.02247 

0.02247 

0.069C& 

0.0690& 

0.03072 

0.037M 

0.03752 

0.03C72 

0.02890 

0.0574S 

0.02371 

0.04021 

0.02890 

0.0574S 

0.02371 

0.04021 

0.01616 

0.02357 

0.01616 

0.02357 

0.01536 

0.03452 

0.01536 

0.03452 


-0.3834098 

-0.2394708 

-0.1558476 

-0.3548763 

-0.0758723 

-0.0276712 

-0.2196882 

-0.1498967 

-0.1624646 


0.004284 


0.05022 

0.03021 

0.0552S 

0.01822 

0.02054 

0.04687 

0.109£" 


42 

42 

q.  7® 


TABLE  4*  (continued) 


43 


i 


m 

i, 


TABIE  4 


o 


Xa 

1.0 

mam 

Xb 

1.0 

0.5 

ex’.. 

0 

-0.0948797 

<*12 

0 

*0.0244703 

0 

-0.0483141 

^21 

°<22 

0 

-0.0497286 

0 

-0.0026621 

Os 

0 

-0.0075401 

31 

0 

-0.0133305 

?32 

«33 

0 

-0.0238455 

0 

-0.0128421 

v(i»i) 

0.001750 

0.002270 

w(0,|) 

0.0005926 

0.0009287 

w(i,0) 

0.001157 

0.001739 

-*5(o,oj 

0.01852 

0.01482 

0.02277 

0.01330 

0.037Q4 

0.06010 

-My(0,0) 

0.02963 

0.050 22 

0.01852 

0.02702 

0.02963 

0.0241^ 

0.03704 

0.03712 

My(0>‘5) 

0.01482 

0.0224^ 

MT 

0.01852 

0.02610 

-Ml 

0.03704 

0.0522g 

*5 

0.01852 

0.02345 

0.03704 

0.03922 

M* 

0.01482 

0.02034 

-MJ 

0.02963 

0.04015 

m*7 

0.01482 

0.01504 

-MJ 

0.02963 

0.02621 

^ = \ 

0.01111 

0.01007 

0.01111 

0.01723 

-*i  = 4 

0.01389 

0.01636 

*6  *8 

0.01389 

0.02065 

M (*) 

0.01852 

0.01351 

0.03704 

O.O3OO5 

<(i5 

0.01185 

0.008975 

-*?(«») 

0.02370 

0.02015 

-0.2337528 

-0.0844420 

-0.1499828 

-0.18910(6 

-0.0155721 

-O.Q357348 

-0.0673447 

-0.1112755 

-0.0706C33 


0.002973 

0.001405 

0.0025pg 

0.02792 

0.01295 

0.09951 

0.08725 

0.03722 

0.01751 

0.Q3716 

0.03222 

0.03520 

0.075CO 

0.02941 

0.04121 

0.02742 

0.056U 

0.01556 

0.02244 

0.008982 
0.0255 2 
0.01942 
0.02912 

0.007452 

0.01922 

0.005155 

0.014C& 


O.OQ4052 

0.002185 

0.003654 

0.03448 

0.01216 

Cx3 

cx5 

0.04919 

0.01260 

0.03954 

0.04456 

0.04635 

0.1202 

0.Q3699 

0.04643 

0.03686 


0.01648 

0.01838 

0.007876 
0.03740 
0.02340 
0. 04037 


45. 


TABU  4.  (continued) 


c 

0.8 

Unite 

A. 

5.0 

2.0 

1.0 

0.89443 

Ab 

6.25 

2.5 

1.25 

1.118Q3 

% 

o<2i 

**22 

«I3 

o<3i 

exit 

♦0.1284045 

♦0.0787999 

♦0.0131284 

0 

♦0.0163288 

+0.0118853 

♦0.0024315 

0 

♦0.Q433052 

♦0.0295710 

♦0.0056153 

0 

♦0.0225577 

♦0.0191622 

♦0.0Q44660 

0 

♦0.0002773 

+0.0004007 

♦0.0001412 

♦0.0006Q40 

0 

♦0.0025828 

♦0.0023419 

0 

♦0.0012468 

♦0.0023937 

♦0.0006341 

0 

cx£ 

♦0.0032919 

♦0.0051678 

♦0.0016686 

0 

CX?2 

33 

-0.0016361 

♦0.0008683 

♦0.0006478 

0 

0.001006 

0.001324 

0.001717 

0.001793 

v(0,|j 

0.0001236 

0.0002894 

0.0005169 

0.0005631 

w(i,o) 

0.0003165 

0.0006902 

0.001144 

0.001230 

>C(i,i) 

0.01211 

0.01534 

0.01901 

0.01967 

*•2 

q*2 

2*2 

qa 

Mj(t,tj 

-Mg (0,0) 

-«;(ofo) 

0.01630 

0.01527 

0.01424 

0.01406 

0.009302 

0.02090 

0.03624 

0.03935 

0.005612 

0.01364 

0.02558 

0.02815 

HS(i.°) 

-Mj(ifO) 

0.005242 

0.01129 

0.01837 

0.01967 

2 

q«2 

0.03840 

0.03396 

0.02903 

0.02815 

vu 

0.03722 

0.03842 

0.03926 

0.03935 

q«2 

Mj(0,*) 

0.003226 

0.007430 

0.01298 

0.01406 

q« 

MT 

0.006731 

0.01216 

0.01851 

0.01967 

q*2 

-Ml 

0.01662“ 

0.02572 

0.03714 

0.03935 

q«2 

m; 

O.Ollcg 

0.01471 

0.01891 

0.01967 

q*2 

“m; 

0.03374 

0.03631 

0.Q3893 

0.03935 

q«2 

m; 

0.006502 

0.009382 

0.01330 

0.01406 

q4 

q«2 

-M2 

0.01548 

0.02012 

0.02678 

0.02815 

Ml 

0.01456 

0.01423 

0.01408 

0.01406 

-Ml 

0.03501 

0.03198 

0.02872 

0.02815 

qa2 

| = | 

0.01279 

0.01179 

0.01074 

0.01056 

q.2 

0.002344 

0.005466 

0.009705 

0.01056 

4 

q 

0.01020 

0.01217 

0.01436 

0.01475 

qa 

-4  = ? 

0.003836 

0.006338 

0.01374 

0.01475 

M(i) 

0.02622 

0.02259 

0.01837 

0.01760 

qa- 

-H*(oj 

O.Q465J 

0.Q4180 

0.03624 

0. Q3519 

qa? 

0.01613 

0.0148? 

0.01298 

0.01259 

q.| 

q« 

"V05 

0.0280Q 

0.02722 

0.02558 

0.02518 

46. 


TiBIE  4*  (continued) 


c 

0.8 

0.5 

Unite 

0:5 

00 

5.0 

2.0 

*b 

0.625 

00 

5.0 

2.0 

«u 

o<i2 

°*21 

*22 

S33 

231 

223 

*32 

*33 

-0.0783882 

-0.0187947 

-0.0396280 

-0.0385983 

-0.0018237 

-0.0060379 

-0.0097162 

-0.0182439 

-0.0093493 

♦0.0394424 

♦0.0017901 

♦0.0336760 

♦0.0073822 

♦0.0000087 

♦0.0062732 

♦0.0002331 

♦0.CQ35093 

♦0.0000093 

♦O.O3O2503 

♦0.0018314 

♦0.0267743 

♦0.0070560 

♦0.00002S4 

♦0.0053878 

♦0.0003675 

♦0.0040001 

♦0.00Q3240 

♦0.0176893 

♦0.0014293 

♦0.0163853 

♦0.0052186 

♦0.0000425 

♦0.0036287 

♦0.0003871 

♦0.0033857 

♦0.0004643 

v(0,|) 

w(i,  o) 

0.002226 

0.0008402 

0.001715 

0.0001584 

0 

0 

0.0003756 

0.0000243 

0.0002309 

0.0006412 

0.0000588 

0.0005133 

1*4/“ 

q*4^ 
q » /to 

-M/(0,0) 

0.02326 

0.00331 

0.05853 

0.04480 

U 

O.OQ4452 

0.009161 

0.0070C2 

0.002832 

0.0086$! 

0.006152 

0.01591 

0.007132 

2 

q«2 

qa2 

S2 

:2*i| 

0.0267S 

0.02352 

0.03942 

0.02044 

• 

0 

0.02072 

0.Q1420 

0 

cJ 

0.008292 

oioz^l 

0.003851 

q*2 

5» 

m; 

-Hi 

m; 

“Hj 

0.026Q3 

0.05220 

0.02382 

0.04131 

Ixd""  . 

0.0Q382S 

0.009822 

0.004341 

O.Olig^ 

0.00837$ 

0.0176S 

0.008612 

0.02122 

% 

q« 

-M< 

Hr 

-HJ 

0.01866 

0.03662 

0.01423 

0.02523 

0.00337$ 

0.008345 

0.009174 

0.0195$ 

0.004174 

0.009500 

0.006549 

0.01800 

0.005314 

0.01122 

0.007818 

0.01614 

q*2 

% 

q« 

0.009668 

0.01562 

0.01685 

0.02041 

0.0075SO 

0 

0.001472 

0 

0.006922 

0.001179 

0.C03894 

0.002788 

0.006152 

0.002838 

0.006812 

0.006184 

S2 

q» 

H (i) 

0.01339 

0.02926 

0.01022 

0.02240 

0.02052 

0.03772 

0.0Q4142 

0.006881 

0.01872 

0.03504 

0.004055 

0.007CS2 

0.01659 

0.03182 

0.0Q3851 

0.007132 

q*3 

qa3 

qa3 

q» 

47. 

T4BIE  4*  (continued) 


I 


♦H- 


48. 


TAB1E  4.  (continued) 


Units 


«<*.*) 

*?(*.*) 

-M5(0f0) 

-*G(0,0) 


-Kj(i,o) 
H*r(0,i) 


i = 5? 


♦0.023^609 

♦0.0015529 

♦0.0219426 

♦0.0061546 

+0.0000311 

♦0.0048069 

♦0.00Q3744 

♦0.0Q38573 

+0.0004169 


0.0006385 

O.OOOQ339 

0.0005108 


0.002237 


0.008386 

0.01792 

0.006716 

0.02342" 

0.004245 

0.009424 

0.007639 

o.oi59£ 

0.0061^2 
0. 0016a 

0.007007 

0.006162 

0.01656 
0.  0313*3 
0.004475 
0.008122 


♦0,0099808 

♦0.0008326 

+0.0097622 

♦0.0032094 

♦0.0000276 

+0.0024Q38 

♦0.0002569 

♦0.0022732 

♦O.OOQ34X8 


0.0009686 

0.00006& 

0.0008611 

0.01406 

0.0068?0 

0.02703" 

0.008560 

0.01385 

0.01388 

0.03092 

O.OQ4444 

0.01390 

0.02802 

0.014CC 

0.03046 

0.005277 

0,01 


0.01361 

0.005152 

0.0Q3304 

0.01070 

0.01036 

0.01385 

0.02703 

0.004444 

0.008560 


0.001174 

0.0000953 

0.001079 

0.01726 

0.006102 

0.Q3452 

0.01220 

0.01726 

0,01220 

O.Q3452 

0.006102 

0.01726 

0.03452 

0.01726 

0.03452 

0.006102 

0.01220 

0.006102 

0.01220 

0.0Q4576 

0.004576 

0.01294 

0.01294 

0.01220 

0.02441 

0.004315 

0.008630 


-0.0122056 

-0.0013975 

-0.0129257 

-0,0051911 

-0.0000766 

-O.OQ37317 

-0.0005528 

-0.0041327 

-0.0007962 


0.001397 

0.0001292 

0.001315 

0.02067 

0.005312 

0.0431,5 

0.017(2 

0.02090 

0.01041 

0.Q38G5 

0.008153 

0.02084 

0.04170 

0.02072 

0.03867 

0.007159 

0.0138g 

0.005521 

0.01074 

0.003970 

0.006175 

0.01533 

0.01573 

0.01045 

O.O2157 

0.004076 

0.008522 


qa4^ 

qa?/D 

qayD 


49* 


TABLE  5.  COHVERGEBCE  OF  A SEQUENCE  OF  APPROXIMATE 
SOLUTICKS  FOB  VARIOUS  SYMMETRICAL  SQUARE  PARKIS 


A«  = «* 


StS! 


♦0.2722222 

0 


0.0013292 

0 

0.02127 

0 

0 

0.04254 

0.004696 

0.009393 

0.01799 

0.03600 

0.01595 

0 

0.03032 

0.06065 


♦0.2801583 

0 

♦0.0510295 

0 

♦0.0326245 


0.0012645 

0 

0.01756 

0 

0 

0.05116 

0.004062 

0.01343 

0.01499 

0.04468 

0.01467 

0 

0.03213 

0.05428 


■ - 1,*,3 

m = 1, 2,3*4 
n = 1,2,3*4 

♦0.2801638 

♦0.2801659 

0 

0 

♦0.0510287 

♦0.0510045 

0 

0 

♦0.0317670 

♦0.0312770 

♦0.0006797 

♦0.0007522 

0 

0 

-0.0026147 

-0.0046450 

-0.0098057 

-0.0168698 

-0.0003257 

0 

-0.0043477 

-0.0114604 

-p.0095106 

0.0012653 

0.0012653 

0 

0 

0.01761 

0.01762 

0 

0 

0 

0 

0.05128 

0.05139 

0.004066 

0.004107 

0.01341 

0.01328 

0.01506 

0.01505 

0.04493 

0.04492 

0.01468 

0.01467 

0 

0 

0.03210 

0.03209 

0.05416 

0.05423 

50 


TABUS  5 


C * 1.0  X*  * X 


n 

a = 1 

B - 1 

mi 

ji 

°*22 

OW 

? 3 

£23 

©<33 

*4 

£24 

% 

40.1108597 

40.0110454 

40.1159889 

40.0109817 

40.0276114 

-0.0003430 

40.0262706 

0.0020850 

0.0007718 

0.0020579 

0.0007727 

-fi(o,o) 
*5  i,o) 
-*S(o,i) 

0*021010 

0.024698 

0.012549 

0.042020 

0.019534 

0.023527 

0.012664 

0.045242 

1 

0.0014262 

0.028524 

0.019675 

0.039350 

0.014194 

0.029734 

0.018484 

0.042826 

0.015758 

0.0092619 

0.015230 

0.0093290 

0.024698 

0.049396 

0.025328 

0.047054 

40.1160046 

40.0109815 

40.0278949 

-0.0003465 

40.0286896 

40.0014561 

-0.0000181 

40.0054417 

40.0014661 


> = 1*2,3»4 
a = l»2t3f4 


40.1160048 

40.0109815 

40.0278906 

-0.00Q3465 

40.0286240 

40.0014329 

-0.0000178 

40.0051624 

40.0003244 

-0.0000587 

40.0000007 

-O.OOQ4054 

-0.0019356 

-0.0022095 


0.0020581 

0.0007726 


0.0020581 

0.0007726 


0.019533 

0.023451 

0.012647 

O.Q45340 

0f 014224 
0.029928 
0.018521 
0.042936 

0.015229 

0.0093315 

0.025294 

0.046902 


0.019538 

0.023455 

0.012647 

0.045344 

0.014224 

0.029921 

0.018522 

0.042926 

0.015228 

0.0093316 

0.025293 

0.046910 


52. 


TABUS  5*  (contintMd) 


c * l«0  Ab  - Xu  = 0*5 


mm 

» = 1 
n = 1 

■ = 1*2.3 
»s  W 

a = 1,2, 3*4 
n = 1,2, 3,4 

felts 

> 
j?  i 

*22 

7f 

?* 

**44 

-0.1424419 

-0.1548472 

-0.1550907 

-0.1550979 

♦0.0283842 

♦0.0286924 

-0,0531865 

♦0.0013214 

-0.0672048 

♦0,0286985 

-0.0552652 

♦0.0013731 

-0.0809059 

-0.0069943 

♦0,0001738 

-0.0292747 

-0.0217654 

♦0.0286986 

-0,0553585 

♦0.0013754 

-0.0817265 

-0.0073321 

♦0,0001822 

-0.0319844 

-0.0288477 

-0.0007988 

♦0.0000198 

-0.0061798 

-0.0102947 

-0.0056857 

0.0Q32713 

0.0033122 

0,0033124 

0.0033124 

9*//® 

v(0,t) 

0.0019634 

0.0019843 

0.0019851 

0.0019851 

q»yD 

<(M) 

0.020606 

0.022424 

0,022477 

0.022472 

JC0.O 

0.063469 

0.068122 

0.068917 

O.O69037 

9*2 

«:i,o 

O.Q31734 

0.030592 

0.03C737 

0.Q3C722 

0.041212 

0.037787 

0.037428 

0.037589 

9? 

*f 

0.029277 

0.028950 

0.028874 

0.028897 

<**2 

q« 

-ft 

0.058554 

0,058204 

0.057554 

0.057416 

0.022322 

0.023790 

0.023722 

0.023714 

0.044644 

0.040340 

0.040196 

0.040214 

0.015455 

0.016156 

0.016164 

0.016161 

q*2 

*^  = «t 

0.023801 

0.023595 

0.023569 

0.023572 

q.2 

*_(*) 

0.015867 

0.015296 

0.015368 

0.015361 

q*3 

q* 

-ij(o) 

0.031734 

0.034061 

0.034459 

0.034518 

£ 2 

f 

<*23 


M 
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-0.597561 

♦0.059538 


O.OO5403 

0.004160 


0.01968 

0.1331 

0.06656 

0.03976 


0.05626 

0.1125 


FOB  A PLATE  SUPPORTED  ONU  BI  COLUMNS 


c • 1.0 


Ak«  A,* 


0.05416 


0.01491 

O.Q4992 


-0.743593 

♦0.064960 

-0.458432 

♦0.017085 

-0.722347 


0.005741 

0.004274 


0.02669 

0.1966 

0.05343 

0.03414 


0.05054 

0.1286 

0.03280 

0.Q3804 


0.01821 

0.04848 


0.005781 

0.004337 


0.02736 

0.2412 

0.05929 

0.02460 


0.05100 

0.1282 

0.03234 

0.03845 


0.01833 

O.Q4725 


* = Exact  Sol'n  Units 

n - l»2v3»4 


-0.7702C79 

♦0.0659718 

-0.5691109 

♦0.0212095 

-1.1627929 

-0.2235778 

♦0.0083322 

-0.7367847 

-0.6860638 


-0.7781116 
♦0.0662663 
-0.6035976 
♦0.0224948 
-1.3126828 
-0.3016386 
♦0.0112414 
-1.0697145 
-1 ..3620243 
-0.1247104 
♦0.0046477 
-0.5051813 
-0.8239495 
-0.6025994 


0.005792 

0.004339 


0.02746 

0.2754 

0.05621 

0.03408 


0.05111 

0.1292 

0.03222 

0.03745 


0. 00.823 
0.04793 


0.005800 

0.004350 


0.02758 

00 

0.05733 

0.02975 


0.05m 

0.1287 

0.03223 

0.03793 


0.01821 

0.04779 


q«yb 
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TABLE  6.  (continued) 


c * 0.6  Xb  ■ Ajj  3 0 


■ w 

■ = 1,2,3 
a = 1,2,3 

Enet  Sol 'a 

tkltc 

°*11 

°*12 

°*21 

3 

-0,3635337 

-0.4610980 

-0.2189531 

-0.3651994 

-0.4297204 

-0.4753682 

-0.2671234 

-0.43U212 

-0.6482044 

-0.0259602 

-0.1748941 

-0.1884084 

-0.4717248 

-0.2527159 

w« 

r*i) 

O.OQ3789 

0.004014 

0.004037 

0.004052 

q4*! 

q»7/J> 

w(C 

0.002004 

0.002140 

0.002156 

0.002185 

w(i,o) 

O.OQ3551 

0.003592 

0.003637 

0.003654 

Sit! 

0.02841 

0.03452 

0.Q3459 

0.03448 

qa2 

0.005955 

0.01112 

0.01311 

0.01216 

qi 

via. 

-M*< 

o,o) 

0.1136 

0.1626 

0.1946 

oO 

-My! 

0,0) 

0.1007 

0.1516 

0.1673 

oc 

vr 

i,o) 

0.05682 

0.04557 

0.05074 

0.04919 

9*2 

;i,o 

0.01191 

0.01874 

0.01127 

0.01260 

q^ 

qa2 

°ti) 

0.05682 

0.04178 

0.04142 

O.Q3955 

mj< 

lo,i 

0.05034 

0.04234 

0.04182 

0.04456 

9* 

M? 

0.05054 

0.04553 

0.04634 

0.04635 

q*2 

qt2 

a* 

-M5 

0.1011 

0.1180 

0.1210 

0.1202 

-Si 

0.03279 

0.Q3781 

0.03699 

0.03699 

0.06558 

0.C4862 

0.C4570 

0.04643 

5 

0.04054 

0.08107 

0.03663 

0.08924 

O.Q3589 

0.08423 

0.03686 

0.08828 

qa2 

9*2 

0.01280 

0.01671 

0.01745 

0.01648 

q« 

-*5 

0.02559 

0.01743 

0.02243 

0.01838 

V 

— wn 

— V 

n.rw./A* 

0-008010 

0.007912 

0.007876 

qa? 

= ■? 

O.Q3775 

0.03803 

O.Q3763 

0.Q3740 

q*2 

q*2 

qa 

-■? 

= >d 

0.02131 

0.02334 

0.02345 

0.02340 

_*5 

^6 

=*s 

0.04261 

0.04087 

0.03998 

C.G4Q37 
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TABIE  6.  (continued) 


c * 0.5 


Parana t«r 


Exact  Sol'n  Unita 


-0.0947776 

-0.1393460 

-0.0354567 

-0.1758046 

-0.0910702 

0.002798 

0.002898 

0.0004096 

0.0005050 

0.002851 

0.002865 

0.03821 

0.04205 

-0.0Q3405 

0.001180 

0.09123 

0.1146 

0.05243 

0.08314 

0.04562 

0,04010 

-0.006810 

0.007961 

0.07642 

0.06148 

0.02621 

0.02769 

0.04398 

0.04116 

0.08796 

0.09933 

0.03935 

0.04217 

0.07871 

0.06734 

0.01967 

0.01806 

0.Q3934 

0.03977 

0.001161 

0.002774 

0.002322 

0.001895 

-O.Of 

o-om/io 

0.01966^ 

0.02337 

0.02866 

0.02942 

0.03421 

0.03330 

-0.1500302 

-0.0584519 

-O.224U03 

-0.1892505 

-0.0134850 

-0.1219897 

-0.0526255 

-0.1897833 

-0.0691311 


0.002906 

0.0005230 

0.002896 

0.04074 

0.001175 

0.1373 

0.1042 

0.04364 

-0.004205 

0.05465 

0.02788 

0.04240 

0.10594 

0.04094 

0.06C73 

0.01809 

0.03759 

0.002740 

0.004061 


0.002914 

0.00053d 

qaT/J) 

0.002900 

0.04068 

qa2 

0.001061 

q*2 

q*2 

oO 

O.Q4245 

2 

q*2 

0.001084 

0.05576 

q«2 

0.02758 

q* 

0.04217 

q<4 

0.10653 

q«2 

0.04117 

q«2 

0.06013 

vr 

q*2 

0.03881 

0.002755 

q«2 

0.002853 

q» 

0.0006890 

oa? 

0.02320 

q«2 

0.02958 

*4 

0.03300 

qi2 
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TABUS  7.  COMPARISON  QF  SOLUTIONS  FOR  PLATES 
FIXED  ON  ALL  FOUR  BOUNDARIES 


(All  Unit*  q«  ) 


e 

Mon.  or  Av.  Mon. 

S - Function 

■ 

Handy*  s 

Flnita  Diffaranca 

Analysis 

Analysis* 

Solution** 

■*Mj(  (0»  s) 

0.0176 

0.0175 

0.0180 

0.0514 

0.0513 

0.0474 

4& : % 

0,0096 

0-0096 

0.0102 

0.0291 

0.0290 

0.0278 

0.8 

0.0092 

0.0094 

0.0173 

0.0178 

O.Q358 

0-0358 

0.Q3U 

0.0424 

0.0425 

0.0404 

4&  * 

0.0049 

0.0051 

0,0199 

0.0179 

♦ 1C) 

0.0097 

0.01Q3 

0.0250 

0.0249 

0.0244 

0.5 

0.00096 

0*00098 

0.00101 

i— vwvi  fry— ra 

0.0101 

0.0102 

0.0142 

0.0143 

0.0099 

-Xyiho) 

0.0207 

0.0207 

0.0201 

♦ 1L) 

O.OOQ49 

0.0005 

0.00053 

-♦ftc  ♦ m») 

0.00789 

0.00785 

0.00558 

tOC  ♦ 1C) 

0.00628 

0.00637 

0.00670 

♦ •%) 

0.01395 

0.0139 

0.0139 

* Huaarical  valuas  art  thosa  coaputad  by  Siass  and  Bsvnark  (1C) 
trca  coafficiants  tabula  tad  by  Vbjtassak  (14). 

**  Area  solutions  nada  at  tba  University  of  Illinois  by  various  parsons 
working  radar  tba  dlraotlon  of  V.  X.  Btwnark  and  C.  P*  Siass.  Each 
span  was  dlvldad  into  sight  aqual  parts;  a a eh  pansl  consistad  of  64 
racbangular  slaaents. 


£ ? iwi  * 
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TABUS  8.  FIRI3B  DUTEHKBCE  S GLUT  I OHS  FOR 
VARIOUS  SIMCTRICAL  SQUIBB  PAULS* 


in 

IB 

1.0 

KB 

5.0 

06 

Units 

0.0067 

0.0037 

0.0029 

0.0023 

0.0018 

0.0014 

q*f/D 

0.0051 

0.0023 

0.0015 

0.0009 

0.0004 

0 

q»Vl 

0.0233 

0.0215 

0.0201 

0.0190 

0.0180 

<4 

0.1897 

0.0665 

0.0410 

0.0233 

0.0102 

0 

9*2 

0.0602 

0.0318 

0.0215 

0.0130 

0.0059 

0 

O.Q335 

O.Q383 

0.0410 

O.Q434 

0.0455 

0.0474 

0.0526 

0.0297 

0.0215 

0.0148 

0.0094 

0.0047 

0.1174 

0.0563 

0.0410 

0.0297 

0.0209 

0.0139 

9*2 

0.0332 

0.0244 

0.0215 

0.0191 

0.0172 

0.0156 

9*2 

0.C430 

0.0409 

0.0410 

0.0412 

0.0414 

0.0417 

q*2 

0 

0.0159 

0.0215 

0.0260 

0.0297 

q*3 

0 

0.0332 

0.0410 

0.0467 

0.0510 

q*3 

* Aron  solutions  aada  at  tho  Tkdvaraity  of  Illinois  by  various 
parsons  working  nadar  tha  dimction  of  V*  M.  Bavnark  and  C.  P#  Siass. 
Each  spaa  was  dlvidad  into  aight  aqual  parts;  aaoh  panal  oonsistad  of 
64  squara  alanants. 
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FIG.  12  POISSON’S  RATIO  CORRECTION  FIG.  13  POISSON'S  RATIO  CORRECTION 

TO  AVERAGE  POSITIVE  MOMENT  TO  AVERAGE  NEGATIVE  MOMENT 

FOR  A SQUARE  PANEL  FOR  A SQUARE  PANEL 


AVERAGE  POSITIVE  MOMENT  FIG.  32  AVERAGE  NEGATIVE  MOMENT 

IE  SHORT  COLUMN  STRIP  IN  THE  SHORT  COLUMN  STRIP 
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FIG.  36  POISSON'S  RATIO  CORRECTION  TO  AVERAGE  FIG.  38  POISSON'S  RATIO  CORRECTION  TO  AVERAGE 

NEGATIVE  MOMENT  IN  THE  LONG  SPAN  NEGATIVE  MOMENT  IN  THE  SHORT  SPAN 


